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Abstract 

A simple lumped parameter based inverse design tool has been developed which provides flow path geometry and 
entrainment estimates subject to operational, acoustic and design constraints. These constraints are manifested 
through specification of primary mass flow rate or ejector thrust, fully-mixed exit velocity and static pressure 
matching. Fundamentally, integral forms of the conservation equations coupled with the specified design constraints 
are combined to yield an easily invertible linear system in terms of the flow path cross-sectional areas. Entrainment 
is computed by back substitution. Initial comparison with experimental and analogous one-dimensional methods show 
good agreement. Thus, this simple inverse design code provides an analytically based, preliminary design tool with 
direct application to High Speed Civil Transport (HSCT) design studies. 

Nomenclature 


A= cross-sectional area 
c= speed of sound 
d= nozzle diameter 
F= thrust 

K= empirical constant 
m= mass flow rate 
M= Mach number 

P,p= Total pressure, static pressure 
R= gas constant 
T= temperature 
S= source term 
u,V= velocity 
Greek Symbols 

p= density 

y= specific heat ratio 
e= small parameter 
X= acceleration potential 
Subscripts. 

1,2,3= primary, secondary, exit (downstream) locations, 
respectively. 

e= exit location, "e" equivalent to "3". 

~= ambient conditions. 

0= total conditions. 



Introduction 


Due to the challenges associated with the High Speed Civil 
Transport Program and the attendant noise suppression issues, 
preliminary design of ejector /mixer nozzles is currently of 
considerable interest. This interest concerning ejector nozzles 
and their noise suppression potential stems from the physical 
mechanisms associated with jet engine aircraft operation. In its 
most fundamental form, jet exhaust noise may be described as a high 
speed "free" jet ef fluxing into a quiescent atmosphere. Lighthill 
(1961) described the fluid mechanics and acoustics of this flow in 
a series of extraordinary papers. Summarizing results from his 
work and others, we may write the rather instructive relationship 
for acoustic power, P a , a measure of "noise": 


P a =K 


p 2 d 2 V 8 

P-C . 5 


( 1 . 1 ) 


where : 

p=density of the jet 
d=nozzle diameter 
p„=density of the ambient air 
c„=velocity of sound 

V=velocity of the jet relative to the surrounding air. 
K=empirical "constant". 

Equation (1) leaves very little doubt that the velocity, V, of the 
jet has a strong influence upon noise generation. Thus, one noise 
"suppression" technique involves minimizing the velocity of the 
jet . 

This strategy of noise suppression must be balanced, though, 
by the thrust requirements of a practical flight vehicle. To 
emphasize this constraint, we consider the ideal nozzle thrust 
relationship (ideally expanded) : 
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F=mV 


(1.2) 


with, m=nozzle mass flow rate. Clearly, to reduce jet velocity for 
fixed thrust, we must increase the nozzle mass flow rate, m. This 
requirement is precisely the situation where a mass augmenting 
device, such as an ejector nozzle, provides a sensible choice. 

Figure (1) presents a representative nozzle attached to a 
turbojet engine and Figure (2) depicts a schematic representation 
of this nozzle for analysis. Fundamentally, an ejector nozzle 
system is merely a mixing chamber in which a high speed primary 
(core) is used to entrain fluid from a secondary flow. The fluid 
dynamic mechanisms associated with this entrainment include both 
viscous and "pressure" components. The actual process involved is . 
extraordinarily complex, involving highly compressible, turbulent 
flow phenomena. More specifically, the entrainment is directly 
related to the local shear layer vorticity (see Townsend, 1976) . 
The so called pressure entrainment is manifested through the global 
momentum conservation statements and dominates long, well mixed 
ejectors. Fortunately, though, rather coarse analyses may still 
be used to provide basic information by integrating between "known" 
locations and "avoiding" complex regions. 

The clear practical application of devices such as the ejector 
nozzle make a hierarchy of predictive analysis of considerable 
interest. This hierarchy may range from simple design correlations 
to state of the art Navier-Stokes simulations. Concentrating upon 
the inverse design problem for which geometry is unknown and 
operating constraints are imposed (as opposed to the analysis 
problem for which geometry is specified and the operating 
conditions are unknown) we propose to develop an simplified, 
mathematical analysis based upon first principles. 

The basic premise of any inverse design analysis, is to define 
a set of design constraints, literally desirable operational 
conditions, and then design a system which meets these constraints. 
From the previous paragraphs two constraints upon the operation of 
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the ejector nozzle system are immediately obvious: 

(1) specification of the nozzle exit velocity, V, which 
we will denote from now on as u e , (V=u e ) . 

(2) specification of the nozzle thrust, F e . 

Specification of the ejector thrust, though, may be at times 
somewhat inconvenient for the cycle analyst who chooses an engine 
of a certain size (primary stream mass flow rate) and iteratively 
looks (via optimization methods) to see if the choice has satisfied 
the thrust requirements. Thus, it will be desirable to include a 
variation of the analysis which is mass flow constrained versus 
thrust constrained. 

These two constraints, exit velocity and thrust or core flow 
rate, are not sufficient to permit computation of the associated 
geometry, but reasonable operational design constraints may be 
imposed. The relationship between constraints and unknown 
variables is most readily discussed by considering the governing 
equations describing the ejector flow. This discussion is the 
basis of the next section. To help gain confidence in the validity 
of this methodology, the report will then summarize available 
theoretical and experimental results. It is hoped, that this 
report will demonstrate the use of this simple analytical in the 
area of preliminary design. 


4 


Analysis 


To provide a design tool of reasonable simplicity while 
retaining sufficient physics to adequately model the flow, the 
inverse design problem is described by a relatively simple set of 
algebraic equations . These relationships comprise both 

thermodynamic definitions and conservation statements with the 
simplifying assumption of quasi-one-dimensional (Q-l-d) flow, 
namely, that all quantities are functions of the streamwise 
coordinate only. Further, the streamwise varying conservation 
equations are integrated (in the streamwise direction) to yield the 
algebraic governing equations. Derivation of these relationships 
may be found in any gas dynamics text, for example, Anderson 
(1982) . These statements imply that the "downstream" conditions 
are fully mixed (cross-stream velocity small and streamwise 
velocity constant, at "mix out" ) an assumption asymptotically 
valid only. We remark that this is probably not an overwhelming 
restriction, in that, we would always strive to design an ejector 
that achieves adequate mixing at the exit. 

Although our inverse design approach has several variations 
(the variations are determined by which quantities we assume to be 
specified and which we assume to be unknown) , most of the governing 
relationships are common to the overall formulation. We will 
proceed to list the governing equations common to this problem and 
then the equations required for the variations. The conservation 
equations (with reference to Figure 2) : 
mass: 

m-L+ihj^ ( 1 ) 


momentum : 


5 


mi u i +p 1 A 1 +ih2U 2 +p 2 A 2 + Jp ( A ) dA =m 3 u 3 +p 3 A 3 


( 2 ) 


and energy: 


*iT 01 -ih 3 T 03 


(3) 


where we employ the following definitions for mass flow at the 
subscripted locations: 

*i s Pi u i A i 

*2 = P2 U 2 A 2 (5) 

*3=P 3 u 3 A 3 (6) 


Note, that the formulation of the energy equation has implicitly 
stated that the specific heat is constant. 

The local ideal gas state relationships written at the three 

locations are: 

Pi-PiR^ (7) 

p 2 =p 2 RT 2 (8) 

and: 

p 3 =p 3 RT 3 (9) 


Further, the definitions of total pressure and total temperature at 
locations (1) and (2) : 
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( 10 ) 


Pi 2 


Ss-(1 (11) 

P 2 2 


and the total temperatures: 


I°l»l+lliM 2 

T x 2 1 


( 12 ) 


T 2 2 * 


(13) 


The Mach number relationships : 


M x 2 = 


2 — U 1 


YRT x 


(14) 


Mj 2 ' — 

12 Y RT 2 


(15) 


and 



(16) 


In the preceding paragraphs, we have developed (16) 
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independent relationships which describe an ejector flow. We now 
proceed to describe the variations of the inverse ejector design 
analysis. By far the most instructive way to proceed is to define 
our unknown quantities as compared to what is known or specified. 
We start with the quantities that are known versus those that are 
unknown regardless of what formulation variation we are concerned 
with. The known quantities include: 

Primary stream (1) : 

• Total Pressure , Poi 

• Total Temperature, T 01 
Secondary Stream (2) : 

• Total Pressure » P 02 

• Total Temperature, T 02 

• Mach Number, M 2 (typically choked, M 2 =l.) 

and the downstream quantities: 

• exit static pressure p 3 =p e 

• exit velocity, u 3 =u e 

where we recognize the exit velocity , u e , specification as a 
manifestation of the noise constraint discussed previously. The 
exit pressure specification is a reasonable design requirement, 
namely, to design the nozzle to achieve ideal expansion. 
Specification of the secondary Mach number is justified by our 
desire to demand as much flow as possible through the secondary 
stream and yet still satisfy the downstream pressure constraint 
without the necessity of . a strong terminal shock, thus forcing 
M 2 =l . 0 . 

The unknown quantities include: 

Primary stream: 

• Mach number, M x 

• density, p 2 

• temperature, T 1 

• static pressure, p 2 

• velocity, u x 

• Cross-sectional area, A x 
Secondary Stream: 
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• density, p 2 

• temperature , T 2 

• static pressure, p 2 

• velocity, u 2 

• Cross-sectional Area, A 2 

• Secondary mass flow rate m 2 
and the downstream quantities: 

• density, p 3 

• Mach number, M 3 

• temperature, T 3 

• Cross-sectional Area, A 3 

• Exit mass flow rate m 3 

Summarizing, we see we have (17) unknowns. Obviously, we have 
more unknown quantities, (17) , than equations, (16) . To proceed, 
we must either reduce the number of unknown quantities or increase 
the number of equations. The choice of what is unknown, depends 
upon what information is available (or conveniently attainable) and 
what quantities we would like to compute. Now, regardless of how 
we constrain the ejector, we must specify the upstream "flow and 
thermodynamic" quantities (such as velocity, Mach number, pressure 
etc.). This is performed in two possible ways: 

(1) Primary Mach number, M x , is specified, bringing our 
equation versus unknown tally to (16) and (16) . This is not 
an easy to obtain quantity and, thus, not normally employed. 

or, alternatively: 

(2) Static pressure matching is demanded between the primary 
and secondary streams. This requirement is desirable to 
reduce "shock" generation due to pressure field imbalances. 
This design constraint takes the form of the simple 
relationship: 


Pi=P 2 


(17) 


Bringing our equation versus unknown tally to (17) versus 
(17). This is the preferable procedure. 
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Now since the number of equations versus unknowns is equal, to 
demand any more information, we must add both a constraint 
(equation) and an unknown for a unique solution to exist. There 
are two computational problems of interest: 

(1) Specification of the ejector net thrust, which introduces 
the relationship: 

Fe SF N =ltl 3 U 3-£ lh l + 'iMU. ( 18 > 


and a correspondingly unknown primary mass flow rate, m x . 

or, 

(2) Specification of the primary flow rate, m x , and the 
correspondingly, unknown net thrust F„. 

Thus, this rather large set of relationships, (1) through (18), 
provides a consistent set of equations describing the ejector 
problem. Although this system is well posed, a detail that has 
been omitted involves closure for the integral within the momentum 
equation. In Appendix (I) we justify the closure to the integral 
which we merely repeat: 

Jp(A)dA«^[p e +^(p 1 +p 2 )l [Aj-A^AJ (19) 


Notice that all terms in this approximation are immediately 
available within the scope of our integral (control volume) 
equations . 

With specification of this term we may now proceed to solve, 
(1) through (19), for the unknown quantities. Now, clearly, 
specification of the thermodynamic /fluid dynamic, (p,p,T,u,M. . •) , 
quantities in the primary and secondary stream is relatively 
straightforward. Essentially, whenever we have access to the total 
quantities and any single static quantity (p, T, M. ..) we may use 
the "local" definitions of the quantities and state equation to 
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compute all of the static quantities. Thus, employing this 
strategy it is convenient to consider all of the primary and 
secondary stream thermodynamic /fluid dynamic quantities as known. 
This computation is not quite as transparent for problems in which 
the closure, equation (17) , is obtained via static pressure 
matching at the confluence and, therefore, bears slightly more 
discussion. Basically, we combine equations (12), (13) and (17) to 
yield: 


goi 

^02 




T 

r=i 



( 20 ) 


which may be solved for the primary stream Mach number: 

Ml 2 = 2 (d+IllMj 2 ) (^)V- 1] 


( 21 ) 


thus permitting us to compute the primary stream quantities. 

The computation of the geometry dependent unknowns which is 
what we are really after (such as, the cross-sectional areas, and 
the mass flow rates) along with the downstream conditions (p 3 , T 3 , 
M 3 ) must now be performed. We relegate detailed derivation of this 
reduction to Appendix (II) and merely quote the results. Although 
explicit formulas are available for the above quantities it is 
probably most instructive to write the governing system for the 
cross-sectional areas (see any good linear algebra text for 
solution method, for example, Noble and Daniel, 1977) . The two 
systems are, of course, different depending upon the problem. 
Writing the system for the "thrust constrained problem": 
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( 22 ) 


(!-£■>* 

U e 

Mi-Pw-^Uo M2 -P w -%u 0 
ifeH 2 -^u 0 u e 


' 



F. 1 

0 


V 


— + s c 

u e 

Pw Pc 


Aa 

= 

F e + S m 

T p 

Y-l e e J 


A 3. 




where the terms are defined (subscripts suppressed): 

V s pu 


(23) 


M = pu 2 +p 


(24) 


fl * C P T °l“y~’[ T °l 


(25) 


The source terms (S c , S m , S e ) found in the RHS column vector are 
introduced to add extra flexibility, for example, skin friction 
losses, S m ; heat transfer, S e ; or mass bleed, S c . Obviously, these 
terms must be specified functions. Throughout the analysis, we 
will assume all source terms to be zero. Further, P w , is defined: 

P w s ^ [P e + ^ (Pi + P 2 > 1 


A few comments about the structure of the system would be 
appropriate. The most "striking" feature of the above system is 
that it is linear . The inverse ejector problem is linear in terms 
of the cross-sectional areas (at least to our level of 
approximation) which is direct contrast to the analysis problem 
which is inherently non-linear in terms of the thermodynamic /fluid 
dynamic quantities. This fact obviously considerably simplifies 
the solution of the problem to the trivial inversion of a (3x3) 
system. 
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In a similar way, the problem may be posed when the primary 
mass flow is specified ( and the net thrust is to be computed) . 

In this case the governing system {also linear) may be written: 


M-P w - 2H 2 ^ ^Pe + Pw 


il-2-L) 


e Y' 

Ho . 2y Pe 


Y“1 U e 


*2 

A, 


2m 1 2l-ihiU 1 +(p w -p 1 )^| 

20 ! 


m. 




u. 


-- 1 ) 


( 26 ) 


with the "auxiliary" equation: 


m, 

1 


( 27 ) 


The matrix system for the primary mass flow rate constrained 
problem is, of course, trivial to invert. We note, also, that the 
"source term" constants have been eliminated for convenience. 

In summary, then, this section has sought to describe the 
theoretical basis of the inverse design methodology. Rather simple 
algebraic relationships of adequate flexibility combined with a 
series of design constraints and assumptions have been used to 
derive a method with the capability of predicting geometric 
parameters {cross-sectional area and mass flow rates) at several 
discrete locations within the ejector. In the next section, we 
will explore the validity of this analysis by comparing to both 
analogous theoretical problems and experimental data. 
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Results and Discussion 


The previous section has sought to outline an elementary 
mathematical analysis describing an inverse design methodology. 
Before this methodology can be applied in any practical manner, the 
range of applicability must be verified. This verification process 
may be conveniently divided into two categories: (1) internal 

consistency checks (such as, simple mathematical constraints and 
cross referencing between analogous, inverse design and simple 
analysis methods); and (2) comparison to available experimental 
data (which is critical since these comparisons provide the only 
available independent verification of the method) . We will begin 
with, the internal consistency of the method since it is 
fundamental and leaves no room for interpretation. 

Due to the rather simple mathematical structure of the problem 
at hand several consistency checks are immediately available to us. 
These might include: 

(1) Conservation statements: the quantities mass, 
momentum and energy are readily shown to be properly 
conserved. 

(2) Mathematical inversion consistency: to guard against 

algebraic "blunders" in the matrix inversions, matrix 
residuals : 

R=[A]x-B (28) 


axe printed and shown to be virtually zero. Although, this is 
certainly a minimal requirement it is a useful test for code 
development /modification procedures . 

(3) "Cross referencing" between variations of the inverse 
design methodology provides another comparison. Since the 
inverse design code analysis has formulation variation that 
are either "thrust constrained" or "mass flow rate 
constrained" opportunities for comparison immediately exist. 
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By supplying these two code variations information from its 
counterpart we have shown these two codes to be self 
consistent in their formulations, 
and, finally, a more stringent test: 

(4) Comparison was made to a forward or analysis integral 
method (see Appendix (III) for a summary description of this 
analysis) . A forward or analysis method may be distinguished 
from the inverse methods described here in that for an 
analysis method, the cross-sectional areas are specified and 
not computed. This comparison was performed for a simple 
constant area case and exhibited adequate correlation. 

As noted, the above consistency requirements, although useful 
and mandatory, are really no more than reflections of the self 
consistency of the mathematics and algebraic manipulation. They do 
not provide any independent information concerning the validity of 
the analysis. To obtain this information, we need to compare to 
either experimental information which, of course, has the greatest 
potential for providing realistic information or comparison to 
higher order numerical analyses (CFD) . 

Limited experimental information is available for a practical 
HSCT ejector nozzle, which is essentially a 2-d, lobed or forced 
ejector nozzle. The model was run in the "thrust constrained mode" 
with thrust set at 52631 lbf (free stream velocity assumed equal to 
zero) . It is probably worth noting, that this rather odd thrust 
value, 52631 lbf, is a manifestation of the iterative modeling 
approach applied by the standard cycle analysis codes. In this 
approach, the thrust is not specified, but iteratively optimized, 
thus, yielding these close to design, 52600 lbf, (but not exact) 
thrust values. Thermodynamic information includes the ideally 
expanded jet velocity (isentropic expansion to ambient pressure), 
Vj=3125 ft/s. The design exit velocity was specified at V mix =1450 
ft/s. This velocity is chosen on the basis of meeting FAR36 Stage 
III noise requirements estimated using semi-empirical noise 
estimation methods (similar to Lighthill's relationship). This 
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relationship is shown graphically in figure (3) . Referring to 
figure (4) excellent correlation is shown between the inverse 
design code and the single experimental mass entrainment value. 
This excellent comparison should be tempered by the fact, that only 
a single experimental value was available which limits our ability 
to more completely verify the method. Further, sensitivity studies 
indicate a strong dependence upon the exit velocity, V,,^, 
specification. Unfortunately, direct measurements of average exit 
jet velocity are not available and we are forced to infer the 
previously stated values. This dependence bears further 
discussion. 

A sensitivity study for a mass flow constrained, generic 
ejector nozzle, using this methodology, is presented in figure (5) . 
These curves exhibit the expected result that more entrainment is 
required to achieve a lower exit velocity. Further, the linear 
relationship between geometry dependent variables, such as, 
entrainment and the thermodynamic relationship is clearly shown. 
We note, that the matrix equations (24) and (28) essentially state 
the same thing. 

Although the previous discussion has centered upon mass 
entrainment (pumping) , the inverse design analysis also has access 
to the required flow path cross-sectional areas at the primary, 
secondary and exit locations. These cross-sectional areas are 
provided for the nozzle for various primary ideal velocities. 
Referring to figure (6), several trends are immediately apparent. 
First, although the size of the primary and secondary flow paths 
change considerably, the overall ejector cross-sectional area 
remains virtually constant. Secondly, streamwise cross-sectional 
area variation is apparently small (relative to length) . This 
confirms that the approximately "constant" cross-sectional area 
ejector is probably a viable design. Both of these trends 
indicate, that although the Mach number of the primary stream is 
relatively high, on the order of Mach 2.0, the ejector geometry is 
dominated by conservation effects. As an example, mass 
conservation is not strongly (in a relative sense) dependent upon 
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Mach number. Consider for example the classical Mach number area 
relationship which is merely a mass conservation statement: 


1 

2 1 

fl + y_1 M 2 ) 

yl 

(29) 

M 2 

Y + 1 

2 / 




Consulting any gas dynamics text, for example, Anderson, 1982 we 
would see that the area ratio, A/A*, varies no more than 1.0 to 2.0 
while the Mach number ranges from 0.3 to 2.25. This insensitivity 
to Mach number and, thus, compressibility causes the cross- 
sectional area trends seen in the previous example. Mathematically 
this means: 


J<J)dA 


♦J 


dA 


pu 

4>=pu 2 +p 
puH 


(30) 


Which is the set of mathematical conditions needed for the 
approximately constant area ejector. 

In any case, estimation of the flow path cross-sectional areas 
represents the first step in the prediction of the ejector 
geometry, and further, an estimate of the cross-section of the 
nacelle itself. Obviously, the actual nacelle cross-section is a 
function of much larger (and more complex) design requirements. 
Work, though, is underway to address these design requirements in 
a rational manner (at a level of fidelity consistent with our 
preliminary design philosophy) . Along these lines, an approximate 
analysis is underway which estimates the streamwise length scale, 
a quantity obviously dominated by the mixing process ( we note that 
our integral formulation has drawn its boundaries or control 
surfaces around this complex region and, thus, avoids such 
complexity at the cost of no streamwise length scale information) . 

A further validation case is available from the fully subsonic 
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ejector experiment performed by (Gilbert and Hill, 1973) . In this 
case a simple, 2 -dimensional, subsonic primary ejector has been 
modeled. Primary mass flow rates are available, therefore, the 
mass flow constrained variation of the analysis is appropriate. 
Necessary, input information includes (for Gilbert and Hill's, run 
number 3) : 

P 01 =4551 . 84 (psfa) 

P 02 =2103.8 (psfa) 

T ol =706 . 0 (R) 

T 02 =533 . 0 (R) 
u e =350 . 0 (ft/s) 
p e =2113.0 (psfa) 

Additionally, the primary and secondary Mach numbers were 
estimated: 

M x =1.0 (choked flow) 

M 2 =.25 (from experimental data) 

Using the inverse design code, estimates of pumping and the exit 
area were obtained and compared to experimental values : 

Given the simplicity of the 
model, these comparisons may be 
considered quite satisfactory. 

To be realistic, it is necessary 
to point out some of the 
uncertainty that enters our 
modeling of this flow. First, 
since complete mixing is 
possible only asymptotically, estimation of the exit velocity is 
difficult. A related difficulty is found in estimating the primary 
and secondary Mach numbers. Note that the experimentally measured 
secondary and primary stream pressure difference: (Pi-p 2 ) /Pi=14.9%, 
is small but not zero. This difference was reproduced in our model 
since we did not demand static pressure matching at the inlet but 
merely imposed the experimentally measured Mach number. Large 
cross-stream pressure gradients would certainly invalidate the one- 



w 2 /w x 

A 3 ft 2 

Exper . 

3.76 

.1278 

Anal . 

4.08 

.1224 

R. Err. 

8.51% 

4.00% 
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dimensionality of our method. In this case, the cross-stream 
pressure gradient appears to have a negligible effect on the 
predicted parameters. Finally, estimation of secondary inlet 
recoveries (primary nozzle recovery was given) may be difficult. 
In spite of these limitations, this example from (Gilbert and Hill 
1973), helps indicate that the inverse design methodology provides 
a useful mathematical model for preliminary design work. 

Finally, since relatively little independent data is available 
to verify the accuracy of this analysis, a comparison to 
equivalent , 1-dimensional, code predictions from other sources is 
desirable. A comparison is shown for a class of generic ejector 
nozzles studied by industry. Although specific information about 
industry methodologies is necessarily incomplete due to the 
proprietary nature of these tools, plausible operating conditions, 
such as, primary and secondary conditions may be inferred. The 
exit velocity constraint, u e , is especially difficult to estimate. 
When the proper exit velocity is applied, comparison is good. This 
comparison is presented in figure (7) for two operating points. 
The comparison is good, but again, this may be a matter of 
consistency between the mathematics of the industry study and our 
method . 

Thus, although data for direct verification is rather sparse 
(hence the need for the analysis) , initial comparisons are 
encouraging. Further, the method has been shown to be both 
internally self consistent- (a minimal requirement) and consistent 
with 1-dimensional predictions from industry. These results are 
probably sufficient to cautiously use the analysis in its intended 
preliminary design role. 


Conclusions 

A simple integral based inverse design model has been 
developed to predict flow path cross-sectional areas and 
entrainment given a physically justifiable complement of design 
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constraints. These constraints include specification of primary 
mass flow rate or ejector net thrust, exit velocity and static 
pressure matching (ideal expansion) . A simple closure to the 
secondary stream static conditions is provided by a primary and 
secondary stream static pressure matching. The resulting systems 
in terms of the flow path cross-sectional areas are shown to be 
linear and, therefore, easily invertible. Comparison to available 
experimental and analogous one-dimensional methods show good 
correlation. Further, the method is easily shown to be internally 
self consistent. Thus, this code provides a preliminary design 
tool for High Speed Civil Transport design issues; and, 
additionally, may provide the basis for more complete design 
methodologies . 

Although this analysis provides an efficient preliminary 
design tool, it is incapable of providing higher information. This 
type of information includes streamwise "mixing length" , pressure 
recovery, flow field information and pumping rates for non-ideal 
mixing. All of these quantities involve the physics of the actual 
mixing process and, thus, may not be modeled by the control volume 
technique employed in this report. A series of higher order models 
based upon ordinary and partial differential equation methods (but 
still much less complex than Navier-Stokes simulations) are under 
development. It is expected that the combination of the simple, 
inverse design model described in this report and these higher 
order models will provide the designer with a simple, yet powerful 
suite of design and analysis codes. 

Acknowledgements 

The author would like to thank Mr. Jon Seidel, Dr. John Lytle 
and Mr. Paul Barnhart for their technical support and input. This 
j'gsgarch was sponsored by the NASA Lewis Research Center, 
Cleveland, Ohio, under contract NAS3-25266. 


20 


References Cited 


(1) Anderson, John D., Modern Compressible Flow, McGraw-Hill, NY, 

1982, 466p. 

(2) Gilbert, G. B. and Philip G. Hill, "Analysis and Testing of 

Two-dimensional Slot Nozzle Ejectors with Variable Area Mixing 
Sections", NASA CR 2257, 1973. 

(3) Lighthill, M. J., Sound Generated Aerodynamical ly, The Baker ian 

Lecture, The Royal Society London, 1961. 

(4) Noble, B., J. W. Daniel, Applied Linear Algebra, Prentice-Hall, 

NY, 1977, 477p. 

(5) Robertson, J. M. , Hydrodynamics in Theory and Application, 

Prentice-Hall, 1965, 652p. 

(6) Skebe, S. A., D. C. McCormick and W. M. Presz, "Parameter 

Effects on Mixer-Ejector Pumping Performance", UTRC Report, 
1988. 

(7) Townsend, A. A., The Structure of Turbulent Shear Flow, 
Cambridge University Press, Cambridge, 1976, 42 9p. 


21 



Appendix I 

Wall Pressure Integral Closure 


The system of governing integral relationships used to 
describe the ejector flow provides a simple, yet physically 
realistic tool for preliminary design computations. This 
simplicity is directly related to the control surface chosen. 
Unfortunately, the momentum equation, since it is conserving a 
vector quantity (linear momentum) , must retain terms that are not 
particularly convenient. This term is the wall pressure or "pdA" 
integral. Closure of the "pdA" integral, in terms of available 
quantities, requires some modeling. 

Instead of merely assuming a convenient form for this term, we 
will briefly introduce a somewhat more sophisticated model, make 
assumptions, and draw what conclusions it permits. First, we 
consider the two dimensional pressure field as described by the 
small disturbance, acceleration potential (Robertson, 1965) : 
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with the attendant "wall" boundary conditions: 
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and, f(x), is the function describing the wall boundaries. 

Now, we may choose to either solve equation (Al.l), (with two 
more boundary conditions, of course) or alternatively simplify it 
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directly. Since our simple integral analysis will be unable to 
"see" two dimensional effects anyhow, we proceed to average by 
introducing the integral operator (essentially integrating (Al.l)) 
in the cross-stream direction) : 




(A1.4) 


Assuming that epsilon is small in equation (A1.2), (which is 
consistent with our inverse design experiences) , simplifying 
equation (A1.4), and recognizing the definition of a cross-stream 
averaged quantity, we may immediately write: 

S.=0 p=c 1 x+c 0 (Al.5) 
dx 2 1 


Note, that this solution has obviously assumed that the flow has 
definite subsonic regions and is therefore, elliptic. 

The purpose of this somewhat lengthy discussion has been to 
"justify" under proper restrictions, the assumption that the 
average streamwise pressure field may be described by the iruear 
function as in (Al.5). With this assumption and the convenient 
specification that the cross-sectional area varies linearly we may 
write the two relationships: 

A(x) =(A e -Ai) (^) +Ai (A1.6) 

and 

p(x) =(p e -Pi) (^) + Pi (A1.7) 

Now digressing briefly, the linear geometry relationship (A1.6) 
bears further discussion for axi-symmetric flow. Considering a 
conical section, (figure 8), we may write the relationship for the 
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radius as a function of the streamwise coordinate "x": 


r (x) =(r 2 -r 1 ) ^+r a — r (x) =^(€^+1) 


€ = 




Thus, 

A(x) =xr 2 [l+c^+€ 2 (-£) 2 ] ■ •Jir 1 2 [l+e^+0(e 2 ) ] 

LLj i-i 


Which justifies, at least to first order, the linear cross 
sectional area relationship, (A1.6). 

Returning to (A1.6) and (A1.7), and eliminating the "x" dependence 
from them we write: 

p (A) = (p e — p^ +Pi (A1.8) 

Thus, substituting into the wall pressure integral (and performing 
the indicated operation) , we may easily write: 

JpdA=^[p e +pJ [A e -Aj (A1.9) 


To complete this analysis, we eliminate p^, by merely assuming: 


Pi 


■5 (Pi + P 2 ) 


(A1.10) 


which is an exact statement when we apply the static pressure 
matching closure assumption. 

Thus, we have obtained a simple algebraic closure 
relationship, which is linear in terms of the cross-sectional 
areas. Although this closure is approximate (viable for small 
divergence angles) , experience shows that this is by no means a 
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severe limitation, since the "solution" cross-sectional area 
variation is small (which might be considered fortunate, in that, 
for large divergence angles, the one-dimensional basis of the 
entire method must fail, due to large cross-stream gradient 
effects) . 
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Appendix IX 

Reduction of Governing Equations 


Previously, governing equations and mathematical results (in 
the form of two simple linear systems) were stated for both "thrust 
constrained" and "mass flow rate constrained" problems. It is 
certainly desirable to fully describe the necessary reduction to 
this point. Due to the formulation, as stated previously, the 
thermodynamic /fluid dynamic quantities (p, u, p, T, M, . . . . ) in the 
primary and secondary stream entrances are essentially known 
quantities. Starting from this "point of view", we quote the 
conservation equations: 

p 1 u 1 A 1 +p 2 u 2 A 2 -p 3 u 3 A 3 = 0 (A2 . 1) 


( Pi^! 2 +p 1 )A 1 +(p 2 u 2 2 +p 2 )A 2 + -^ (3p e -^ (P 3 + P 2 ) ) (A 3 -A 3 -Aj ) =p 3 A 3 + p 3 u 3 2 A 3 

and 

P i u iAiC p T 0 i + p 2 u 2 A 2 c p T 02 -p 3 u 3 A 3 c p T 03 (A2 . 3 ) 

We note, that u 3 =u e and p 3 =p e (for any of the inverse design 
variations) are both specified quantities. 

We now proceed to consider the "thrust constrained" problem. 
For this variation, the required net thrust of the ejector is 
assumed specified: 

Fe =F N = P 3 u 3 2 A 3 -u.(p 1 UiA 1 +p 2 u 2 A 2 ) = const. (A2.4) 


Thus, we may simply eliminate the term on the right hand side of 
momentum : 
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p 3 u 3 2 A 3 =F e +u.(p 1 u 1 A 1 +p 2 u 2 A 2 ) 


(A2.5) 


and, analogously the right hand side of the mass conservation 
relationship becomes : 


P3 u 3 A 3 = 7T + ^T-( P 1 u 1 Ai + P 2 u 2 A 2 ) 

U e u e 


(A2.6) 


Now, reduction of the right hand side of the energy equation 
into either known or unknown cross-sectional area terms, requires 
somewhat more effort. We begin by solving (A2.5) for the exit 
density, p 3 : 


P 3 = 


F e 

A 3 u e 2 


+u.( 


Pl u i A l + P;UA } 
u e 2 A 3 u e 2 A 3 


(A2.7) 


and by state: 

(A2.8) 

3 p 3 r 


and, finally, by the Mach number relationship: 


m 3 2 = — tF e +u.(p 1 u 1 A 1 +p 2 U 2 A 2 ) ] 


(A2.9) 


Thus, we write (with some algebraic simplification) : 

p 3 u 3 A 3 C p T 03 =-I I p eUe A 3+ l Ue F e+ 1 p^u^u^^ P 2 u 2 u - A 2 u e (A2 . 10 ) 


At this point, it is clear, that the above relationships are all in 
terms of either known quantities or the unknown cross-sectional 
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areas. Collecting terms, the linear system (24) is recovered. 

The alternate formulation, "mass flow rate constrained 
problem" is derived analogously. We proceed by computing the exit 
location density via state: 


_ YPe 
P3 yRT 3 


(A2.ll) 


and substituting into the momentum equation: 


A l U lPl A l + Pl U 2 2A 2 + P 2 A 2 =P e A 3 + 


YPe 

YRT 3 


U o 2A 3 


(A2.12) 


The next step is to solve for the speed of sound using the energy 
equation placed in the form: 


*iT 01 + PjUjAjTo* = ^u e A 3 ( 1 + 1^ ^ 


(A2.13) 


to yield: 


YRT 3 = 



^lToi + P2 U 2 A 2To2 ^ ^ 
pm 


(A2.14) 


Back substitution into energy and substitution into momentum, 
(collecting terms) yield the matrix system, (28) . Thus, we have 
arrived at the desired linear systems in terms of the cross- 
sectional areas. These systems represent the "fundamental" 
mathematical result of our inverse design analysis. 
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Appendix III 

Summary Derivation of an "Analysis" Mixing Problem 


As a simple test of consistency the "inverse" formulation was 
compared to the "forward" or analysis formulation. Although the 
analysis presented is restricted to constant area mixing 
(recognizing that this could be easily relaxed using the pressure 
closure relationship derived in Appendix I) , it provides a useful 
test of the inverse design methodology. Further, it also 
illustrates the contrast in solution structure (notably non-linear 
versus linear) exhibited by the analysis method. This appendix 
briefly summarizes this formulation. 

The governing equations ( quasi-one-dimensional , integral 
relationships) are: 

ihj+rt^ =m 3 sW ( A3 . 1 ) 

lh 1 U 1 +ft 2 u 2 +p 1 A 1 +p 2 A 2 =ih 3 u 3 +p 3 A 3 ( A3 . 2 ) 


and 


I * l iT 01 +ri'i 2 T 02 -m 3 T 03 E 


(A3 . 3 ) 


Now, recall that for the analysis mixing problem geometry cross- 
sectional areas and the both upstream quantities are given (hence 
the definition of W,E and P) . Our task in this problem is to 
compute the "fully" mixed conditions (subscript 3). Accordingly, 
we eliminate pressure and density through the state and mass 
conservation relationships : 


29 



(A3. 4) 


P 3 =P3 U 3 R TT =W 

U 3 


1 

(YRT 3 ) t 1 

y a 3 m 3 


and noting: 


m 3 u 3 =W(yRT 3 ) X M 3 


(A3. 5) 


With these quantities specified, we may substitute into the 
momentum equation: 




(A3. 6) 


Y ML 


rewriting and squaring both sides: 


M 2 ( £n;om ) 2 y 2 =YR1 p 3 [ ± +yM3 2 ] 2 

w 


(A3. 7) 


Now, eliminating the speed of sound through the energy equation: 

yR— 

yet,- v " < A3 - 8 > 

1+IJ:M 3 2 


we may write (after collecting terms) the non-linear (4th order 
relationship) for the exit Mach number: 


[ y 2 (1-^)+^]M 3 4 +y(2-G)M 3 2 +1=0 


(A3. 9) 


where "G" is the non-dimensional grouping: 
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(A3. 10) 



REW 


The solution of the "4th" order non-linear equation is really 
quite simple when we recognize that the substitution M 3 2 =Z, 
immediately reduces the equation to a quadratic with two roots . 
These two roots denote the supersonic and subsonic solutions of the 
mixing process (the negative roots are obviously trivial) . It 
would go well beyond the scope of this Appendix to discuss this 
solution further. The other quantities, velocity, temperature, etc 
are available by back substitution. Although this analysis did not 
play a direct role in the inverse design problem, it did provide a 
useful check, and point of reference for the inverse methodologies. 
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Figure 2. Ejector Schematic. 
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Vmix (ft/sec) 

Figure 3. FAR36 Stage III Requirements. 



Videal (ft/sec) 


Figure 4. Inverse Design Analysis/Experimental 
Comparison. 


33 







Figure 7. Comparison to Analogous, Industry 
Prediction Methods. 



Figure 8. Linear Axi-symmetric (Conical) Geometry. 
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